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\f} ■ I. INTRODUCTION 

<n : 

Among the many problems which attend the canonical quantization of general relativity is the fact that all gauge- 
invariant quantities, which might qualify as observables, necessarily commute with the Hamiltonian, and are therefore 
■ constants of the motion. This intriguing problem has received wide attention, which has been comprehensively 
reviewed, for example, in [jjjpl (see also It seems possible to obtain some insight into the "problem of time" from 

simple quantum-mechanical models which share with general relativity the essential feature of time-reparametrization 
ON ' invariance, and the purpose of the present work is to explore in some detail the interpretation of time evolution in 
one such model. 

It is widely thought that the time evolution of physical quantities which is apparent to actual observers in our 
universe should be understood, in some sense, as evolution relative to the state of a physical clock, which provides 
an observer's local definition of time. A particular implementation of this idea has been proposed by Rovelli 
who defines "evolving constants of the motion" by eliminating the arbitrary time coordinate t in favour of the value 
Jh . of a physical clock variable. Thus, in ||, he considers a model of two harmonic oscillators having the same frequency, 
whose equations of motion can be solved for two functions, say xi(t) and x 2 (t). Neither of these functions gives a 
^ . gauge-invariant description of time evolution, but by eliminating t, Rovelli arrives at a new function X\(x2) which is 
gauge invariant (or almost so - see section 4 below) provided that x 2 is regarded as a real parameter, rather than 
as a dynamical variable. Classically, at least, this defines a family of observables, £i(s), which can be interpreted as 
the value of x\ when x 2 — s. However, the quantum-mechanical operator corresponding to x\(s) is approximately 
self-adjoint (and approximately obeys a standard Heisenberg equation of motion) only when restricted to a subspace 
of the physical Hilbert space which Rovelli calls the "Schrodinger regime" . Since this is a large-quantum-number 
regime, he suggests that the notion of time evolution does not exist at a fundamental level, and emerges only as a 
semiclassical property of macroscopic objects. 

In this paper, we study a model whose dynamical variables are an oscillator and a free particle, the latter serving 
as a clock whose reading is (classically) a measure of geometrical proper time. The classical and quantum dynamics 
of this model are discussed in Section II where, in order to obtain a well-defined quantum theory, the free particle 
is treated as the low-frequency limit of a second oscillator. We find it possible to define a parametrized family of 
observables x(t) evolving with a parameter r which classically coincides with the geometrical proper time. The 
linearity of the clock variable turns out to be inessential for this purpose. Moreover, this evolution is governed by 
an exact Heisenberg equation of motion, and is not restricted to a macroscopic regime. The introduction of a clock 
variable permits time evolution by relaxing the energy constraint on the oscillator, and indeed the energy of the clock 
can saturate this constraint. In Section III, we investigate the possibility of a "test clock" by restricting the state of 
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the clock to a narrow range, which approximately reproduces the constraint on the oscillator alone. This restriction 
imposes a limit on the resolution with which the value of x(t) can be determined, but there is a semiclassical limit in 
which both the energy of the clock and the eigenvalues of x(t) are fairly sharply defined. The interpretation of time 
is discussed in detail in Section IV where, in particular, we point out that the clock cannot itself be regarded as a 
physically observed object, but should rather be taken as representing the observer's intrinsic sense of time. 



II. CLASSICAL AND QUANTUM DYNAMICS OF AN OSCILLATOR AND CLOCK 

Consider the Lagrangian 

L = A (x 2 + q 2 ) ~ {u 2 x 2 - 2E ) , (1) 

defined on a (0 + l)-dimensional manifold with "time" coordinate t. The only geometrical variable is lapse function 
A(i), which we take to be strictly positive. The action J Ldt is easily seen to be invariant under a change of time 
coordinate, with x'(t') = x(t), q'(t') — q(t) and A'(i') = (dt/dt')X(t), where t' is any increasing function of t. The 
equations of motion 



Id fldx 
Xdt \X~dt 
ld_ ridq 
Xdi \X~di 

obtained by varying x and q have the general solution 



-u?x (2) 

(3) 



x(t) — oi cos-0(i) + <i2 sin ip(t) (4) 

q(t) = a + c 2 i/j(t) , (5) 

so that the clock variable q is linear in the proper time w -1 ?/>(£) = f Q X(t')dt' . The constraint 

^ 2 1 ^ +u 2 x 2 = 2E (6) 



obtained by varying A implies the relation 

C = a\ + al + cl = pl , (7) 

where p\ = 2E /lo 2 . 

For this model, the solution space (ignoring the factor A of gauge functions A(t), all of whose points are gauge 
equivalent) is S 2 x R. With the parametrization 

ci\ = po cos x cos 7 (8) 
«2 = -pocosxsin7 (9) 
c 2 = posinx , (10) 

the angles \ and 7 (— ir/2 < x < 7r/2, < 7 < 2tt) provide coordinates on S 2 , while c\ is the coordinate on R. The 
general solutions (|J) and (||) can be written as 

x(t) = p o cosxcos(-0(i) +7) (11) 
q(t) = ci +p sm X i>(t) (12) 

and the presymplectic form on the solution space is 

fi = uj(dai A da,2 + dc\ A dc 2 ) 

= \ u Po sin ( 2 x) d X A d7 - ujp cos x^x A dci . (13) 
As can be seen from ([ll]) and ( |l2| ) , time evolution is generated by the Hamiltonian vector field 
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X c = |- + p sin X A (14) 

associated with the constraint C. As expected, we find that Xq annihilates f2, ix c ^ = 0, so the reduced phase space 
is found by identifying all points along each gauge orbit, which is an integral curve of Xq. Equivalently, points of the 
reduced phase space correspond to distinct gauge orbits. The angle x is constant along each gauge orbit. For each 
value of x (except x = 0, ±7r/2), the orbits are helical curves on the cylinder S 1 x R coordinatized by 7 and ci, with 
a period along R of 2irpo sinx- Distinct orbits can therefore be labelled by the values of c\ in an interval of length 
27rposinx, the end points of which are identified, since they belong to the same orbit. The set of distinct orbits is 
thus topologically S 1 . For x = 0, the gauge orbits are circles on the cylinder, and there is a distinct orbit for each 
ci e R. For x = ±tt/2, the amplitude of oscillation in ([ll]) vanishes, so all values of 7 are equivalent. In this case, 
the S* 1 of gauge orbits collapses to a point. 

We see that the reduced phase space is the union of three sets: 

(i) the open disk D + , which is the direct product of the interval < x < tt/2 with the circle S 1 , the circle at 
X = 7r/2 being contracted to a point. Physically, states in this region are those in which the clock runs forwards. 
The gauge-invariant observables which distinguish the states are Xj which sets the amplitude of the x oscillation 
and also (via the constraint) the rate of the clock, and ci, which characterises the phase of oscillation at which 
the clock reads zero. 

(ii) the open disk _D„, which is similar to D + , but corresponds to the interval — ir/2 < \ < and contains those 
states in which the clock runs backwards. 

(hi) the line Rq, with c\ £ (—00, +00) and x — 0. In this region, the oscillation has its maximum amplitude, while 
the clock has a constant reading ci, which is the gauge invariant observable distinguishing the states. 

This reduced phase space is not a manifold - a fact which presents serious difficulties for the method of geometrical 
quantization, and we have not succeeded in obtaining a fully consistent quantum version of this model. A quantum 
theory of a related model will be described below. Within the classical theory, we can construct gauge-invariant 
observables which evolve with proper time, as measured by the clock, by restricting attention to the subspace D+ UD_ . 
Consider that region of the solution space where \ 7^ and identify the points (x,7,ci) and (X:7>ci + 27rposinx). 
On the resulting space, we replace the coordinates c\ and 7 with 

c=^— (15) 
Po sinx 

7 — — ) mod 2-7T . (16) 

Po sinx/ 

In terms of these coordinates, the generator (H) of gauge orbits is 



and the presymplectic form (|13[) becomes 



n=-up 2 sm(2 X )dxAda. (18) 



As expected, this Q is independent of c, and is also the symplectic form on the region D + U £>_ of the reduced phase 
space, now coordinatized by x and a. 

With the coordinates (x, a, c), the solutions to the equations of motion read 

x(t) = po cos x cos (ip(t) + a + c) (19) 
q(t)=p sin X ty»(*)+c) . (20) 

The quantity t = (ip(t) + c) /u) can be identified as the proper time which has elapsed since the clock read zero and 
(with an obvious economy of notation) the value of x at this time is given by 

x(t) = po cos x cos(wr + a) = X cos(ojt) + uj~ 1 P sin(wr) , (21) 

where 
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X = po cos x cos a (22) 
P = — uipo cos x sin a . (23) 

Since they do not depend on c, the quantities X and P are gauge invariant. Thus, if r is regarded as a parameter, 
rather than as standing for the expression (ip(t) + c) /u>, then x(r) is a gauge-invariant function of r. Classically, this 
function can be interpreted as "the value of x at proper time r" , given that r = is located at the event q — 0. It is 
similar to the "evolving constants of the motion" defined by Rovelli fq|| but with an important difference: whereas 
Rovelli's evolution is with respect to the actual value of a physical clock variable, r refers (at the classical level) to 
the geometrical proper time. The symplectic form can be expressed in terms of X and P, with the satisfactory result 

£1 = dX A dP , (24) 

so X and P have the canonical Poisson bracket algebra. Moreover, the evolution of x(t) is easily seen to be governed 
by the usual equation of motion 



dx(r) 
dr 

where 



{x(t),H } , (25) 



!/d2 , , .2 V 2N _ 1 ( dx ( T ) 
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H„ = :(P 2 +^) = :M^] l-^.nrrj il>.i) 

is the Hamiltonian for the oscillator alone. 

To obtain a quantum version of this model in a controlled manner, we consider the Lagrangian 

L = Y\ { - ±2 + ^~\ ( wV + N ~ Wq2 ~ 2E ^ ' (27) 

where N is an integer, with a view to recovering the model (Q) in the limit N — > oo. Somewhat surprisingly, the 
results are essentially independent of iV. The classical analysis is similar to that given above. If the solutions to the 
equations of motion are written as 

x(t) = ai cos^i(i) + a-i sin^i(i) (28) 

q(t) = h cos (N- 1 ^)) + b 2 sin (JV - V(*)) , (29) 

then the quantities 



a = ^/ — (ai + ia,2) (30) 

■^Qb-iln) m) 

have Poisson brackets {a, a*} = {fe, b*} = —i, which will shortly be promoted to commutators. The constraint 
a\ + + N~ 2 (bl + b\) = Po can be solved in terms of the coordinates ol and c, replacing (|l9| ) and (|o|) with 

x(t) = po cos x cos (ip(t) + a + c) (32) 
q(t) = Np QS m X sm(N- 1 (^(t) + c)) . (33) 

For any finite value of N, the reduced phase space of this model is S 2 , the angles x an d a (0 < \ < n/2, < a < 2ir) 
providing coordinates such that the points x — 0: 7r /2 are opposite poles. In particular, for x — 0j the angle a in 
( ^2| ) is indistinguishable from c, so all values of a are gauge equivalent. Each point of this phase space, of course, 
corresponds to a periodic oscillation, for which the values ip(t) — 0, 2iV7r may be identified. The limit N — > oo is 
actually singular, in the following sense. To achieve a finite limit in ([33]), with % ^ 0, we may assume that (if>(t) + c) 
either has a finite value, or differs by a finite amount from iV7r. In the limit, therefore, each point of the phase space 
S 2 represents two distinct solutions, corresponding to the half-periods of oscillation in which q is an increasing or a 
decreasing function of ip(t). Clearly, the two limits of (|3^) differ by a sign, which can be regarded as the sign of x, 
and correspond to the two regions D± of the phase space described above. To reproduce the region Rq , we must first 



take N — > oo with x 7^ and then take x ~~ * 0, with c given by (15) 
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For finite values of N, the reduced phase space S 2 is a manifold, and the model can be quantized straightforwardly. 
However, since this manifold includes the point \ = 0> we may anticipate some difficulty in interpreting the variables 



X and P (equations ( 22 ) and fl23)) as operators on the physical Hilbert space. Since all values of a are gauge 
equivalent at x — 0, X and P (and hence also x(t)) are not gauge-invariant functions on the reduced phase space, 
although they are gauge-invariant functions on any region which excludes the point \ — 0- I n anv such region, the 
expressions (Q) and (gj) for the symplectic form are still valid, so X and P still have the canonical Poisson algebra. 

Models of this kind can be quantized in several more or less equivalent ways (see, for example P,[s|-[ic| ) ■ The essential 
result is that one can identify gauge-invariant functions (si, 52,33) = (wpo/4)(sin(2x) sina, cos(2x), sin(2x) cos a), on 
the reduced phase space, whose Poisson algebra, {st, Sj} = CijkSk is the Lie algebra of SU(2), and the Casimir invariant 
• s i + s 2 + s i is proportional to the constraint C. The physical Hilbert space is therefore an irreducible representation 
of SU(2) whose dimension is set (up to operator ordering ambiguities) by the value of E . For our purposes, it is 
convenient to construct this Hilbert space according to the Dirac prescription, taking a and b as the basic variables 
(a similar route is followed in p0|). 

On promoting a and b to quantum ladder operators, with the usual commutators [a, a^] = [b, b^] — 1, [a, b] = 
[a, b^\ = 0, we obtain an unconstrained Hilbert space, spanned by the vectors |m, n), with 



a|m, n) — \/m\m — 1, n) , a' \ra, n) = \Jm + l\m + 1, n) (34) 



b\m, n) = \/n\m, n — 1) , 6 T \m, n) = \/n + l|m, n + 1) . (35) 

The physical Hilbert space is the subspace of vectors satisfying the constraint 

(a^a + N^rfb)^) = v\i/>) , (36) 

where v = lopq/2 = Eq/uj. Here it is assumed that any constant arising from factor ordering in the constraint operator 
has been absorbed into E a . Quantization clearly requires that Nv be an integer, and we define 

v = v + u'/N , (37) 

where v and v' are integers, with < v 1 < N. This physical Hilbert space is spanned by the vectors 

\m)) = \m,N(v-m) + v') ,0<m<v. (38) 

We would now like to realize X and P as gauge-invariant operators, acting in the physical Hilbert space, but, as 
anticipated, this is not quite straightforward. Classically, we can define variables A and A*, with Poisson bracket 
{A,A"} = -i,b, 

A-J%(x + i.p)-*(-£ f Y . (39) 



2 V w J \Vb*b 
Quantum-mechanically, the operator ordering 

N 



A = a[{tfb)- 1/2 b ] } (40) 
Al = at (b^b)- 1 / 2 ) (41) 



ensures that these operators have the expected properties A\m)) — ^Jm\m— 1)) and A^\m)) = y/m + l|m + l)), except 
that the action of A' on the maximal state \v)) is not well defined. To proceed, we introduce the regularized operators 



A e = a 



((6t6 + e )-V2 6 t) (42) 

A\ = a) (bitfb + e)- 1 ' 2 } . (43) 

These are well-defined, gauge invariant operators, and we find in particular that -A||z/)) = 0. We can now define 
operators A and by 

A\m)) = lim A e \m) = y/m\m - 1)) (44) 

e— >0 

A f |m)) = lim A\\m) = Vm + l(l - S mv )\m + 1)) (45) 

e— >0 



5 



These operators are also well defined and gauge invariant, but they have the anomalous commutator 



[A, At] = 1-0, 

where 9 projects onto the maximal state: 9\m)) = (1 + v)8 m ^\m)) . The operators 



X = 
P = 



2cj 



1/2 



(A + At) 



(?f *) 



(46) 

(47) 
(48) 



are well defined, gauge invariant and self-adjoint, but they have the anomalous commutator [X, P] = i[A, A'] = i(l—8). 
The anomalous term is nonzero only when acting on the maximal state \u)), where the clock has the smallest energy 
allowed by the constraint (|36|). Classically, the interpretation of X and P as the position and momentum of the 
oscillator "at time r = 0" is ambiguous in the state x — 0, where the clock permanently reads zero, and therefore 
does not distinguish this instant of time, and the anomaly can be understood as reflecting this fact in the quantum 
theory. Formally, the classical versions of (E2h and (|43|) correspond to regularized variables 



X e = po cos x cos a 



1 



2c 



NujPq sin x. 



-N/2 



(49) 



and P t = —Lupota,naX e . These are truly gauge-invariant quantities, having the unique values X e = P t = at x = 0, 
but do not have the canonical Poisson algebra. However, they differ significantly from the original X and P only 
where sin 2 x is n ot much greater than t/ujp^, and when e is sufficiently small, this is a small neighbourhood of the 
point x — 0. 

It is simple to show that 



9 A = A^9 = 
[A, A f A] = A 
[A f , A f A] = -A f , 

and it follows from the latter relations that the Heisenberg equation of motion 



dr 



i[H ,x(r)] 



(50) 
(51) 
(52) 



(53) 



reproduces ( pl| ) in the expected way, provided that the Hamiltonian has the factor ordering Hq = ljA^A. 

At this point, we have a formalism which, taken at face value, is equivalent to ordinary time-dependent quantum 
mechanics. There is a physical Hilbert space, spanned by the vectors |m)), a set of operators (X, P) on this space and a 
Hamiltonian Hq which generates time evolution through the standard equation of motion ([53]) . Whether this formalism 
should be taken at face value is, of course, another matter, and the interpretation of the model will be discussed in 
Section IV. Apparently, one can ask, and answer, questions such as "given that x was determined to have the value x\ 
at time t\ , what is the probability that it has the value Xi at time T2?" . Acording to standard quantum mechanics, this 
question is legitimate only if x\ and xi are eigenvalues of x(t±) and x(t2). Because the classical range of x is restricted 
by the constraint and, correspondingly, the physical Hilbert space is finite- dimensional, these eigenvalues form a finite, 



discrete spectrum, which is easily found. Using (^ 
properties (E3) and (E5I) then imply that 



we can express x(t) as x(r) = (2ujy 1 / 2 (Ae~ iuT + A ] e iuJT ). The 



\x jt T)) =Nj y2(2 m m[)- 1 / 2 H m ( l j 1 / 2 x j )e imUT \m)) 



(54) 



is an eigenvector of x(t) with eigenvalue Xj, where H m is the Hermite polynomial, Nj is a normalizing constant and, 
in order to truncate the series at m = v, uj^l^Xj must be a zero of H v+ \. Since there are v + 1 of these zeros, the 
states ( p4[ ) span the (v + l)-dimensional Hilbert space. Also, since p(r) = x(t) = lox(t + n/2u>), the eigenvalues of 
p(r) are u>Xj. 

The quantum dynamics of this model can be illustrated by considering the coherent state 



\z)) = C{z)e 



m=0 



(55) 



G 



where C[z) is a normalizing factor. The time-dependent wavefunction for this state is 

•4>(xj,T;z) = {{Xj,t\z)) 

(56) 

where z(r) = e~ liA)T z. If z — y/uj/2(x + (i/uj)p ), then z(t) — ^/w/2(xo(r) + (i / uj)p (t)) , where (xq(t),Po(t)) is the 
classical trajectory passing through (xq,po). For an unconstrained oscillator (corresponding to the limit v = oo), the 
coherent state is, of course, a Gaussian wave packet, with |^| 2 oc exp(— uj(x — -\/2/cj5Rz) 2 ), whose peak follows the 
classical trajectory. The wavefunction ( |56| ) is defined only at the discrete values Xj, and is not simply a function of 
Xj — z(t). Nevertheless, if xq is well within the range defined by the largest and smallest Xj then, even for quite small 
values of v, these discrete values follow the Gaussian packet rather closely, as illustrated in figure 1. 



III. TEST CLOCKS 



In classical general relativity, one can assess the physical characteristics of a spacetime by examining the trajectories 
of "test particles" which follow time-like geodesies, but do not contribute to the stress tensor. Here we consider the 
possibility of a "test clock", which might be used to reveal the time evolution of a quantum universe, without itself 
contributing significantly to the quantum dynamics. Whether this notion is ultimately meaningful or useful depends 
on how the quantum theory is to be interpreted, and this is discussed in the following section. 

Suppose, specifically, that the clock variable q is deleted from the Lagrangian ([!]) or (p7f). Classically, there seems 
to be a sense in which the state of the remaining oscillator evolves with proper time as x(t) = y / 2Eo/uj 2 cos(wt), 
although x(t) cannot be regarded as a gauge-independent observable. In the quantum theory, the Hilbert space 
contains only a single state, from which nothing corresponding to this apparent time evolution can be extracted. By 
adding the clock, we are able to obtain a genuine observable x(t) and a multi-dimensional Hilbert space on which 
proper time evolution can be represented in a gauge-invariant manner. However, the energy constraint which applied 
to the original oscillator is relaxed by the presence of the clock, which can itself saturate the constraint. The idea 
now is to see whether the apparent time-dependence of the clockless model can be revealed by restricting attention 
to those states in which the clock has only small energies, so that the energy constraint of the clockless model is 
approximately realized. It may be anticipated that this restriction will impose a limit on the resolution with which 
values of x(t) can be determined, and we wish to examine the nature of this limit. 

If there is to be a non-trivial time evolution, several states must still be available. We therefore consider a (loosely) 
semiclassical situation, with v very large, taking a band of states with m ~ v to be actually available for use. Our 
calculations will now be approximate, and we begin by finding an approximation for the eigenstates (|54|). It is 
convenient at this point to use the Schrodinger picture, and deal with the eigenstates \xj, 0)) of X — x(0). When v is 
large, the v + 1 zeros of H u+ i corresponding to the eigenvalues of x(t) lie roughly between — \[2~v and \/2z/ and the 
spacing between them is therefore of order z/ -1 / 2 . We will take these eigenvalues to form a continuum. For large m, 
the asymptotic formula for H m Q , 



Hm{z) * ^ rjly+i) cos (( 2m + 1 ) 1/2z " m ^ 2 ) ( 5? ) 



yields 



(2 4n+r (4n + r)!) 1/2 H in+r (z) « {2ir)- 1/4 e z * /2 n - 1/A c r (Vs^z) , r = 0- • • 3 (58) 



where c r (9) = cos(# — tti/2). Taking v to be of the form Av 1 + 3 and m — Av'a + r in (|34j), replacing the sum on a 
by an integral, and defining 



we obtain eigenstates of £ in the form 



10) = iff fda a^'^^^OM) 

" n J° r=0 



(60) 
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which are orthonormal in the limit v — > oo. 

Ladder operators = A and A^ 1 ^ = A^ are realized in this aproximation by 

aM \o) = vj^U-^)\o) (6i) 

AM\a,r)) = ^ (a 1 / 2 - ^ 1/4 

\ v oa 

x\a,r-fi)). (62) 

Because the maximal state annihilated by A' now has zero measure, these operators have the canonical commutator 
[A, A*] = 1. For the same reason, the coherent state is now an eigenstate of A. Corresponding to the rescaled variable 
£ @ we define 

A\z)) = ^z\z)) (63) 

and obtain the usual wavefunction ip(£, z) — ~ exp(— t/(£ — z) 2 ). 

Supposing that the states actually available for use are those in a narrow range of a, we introduce a window function 
/(a) and define approximate eigenstates of £ by 

\i))f = f 1 da a-y*f{a) V c r (2zV^)K r» . (64) 
In the same way, we can define an approximate coherent state 

pl 3 

\z))f = d Q «- 1 /V(«)Vi(a^)|a,r)), (65) 

J ° r=0 

where 

tl> r (a,z) = J d£ip(£,z)c r (a) . (66) 

As in ( |56] ) the time evolution of this state is obtained by making the replacement z — > z(t). 

To obtain approximate analytical results, we define 7 = ^/a and choose a Gaussian window function f("f) — 
/oexp(— (7 — 7 ) 2 /2A), where /o is the appropriate normalization factor. We choose 70 to be a little less than the 
maximal value 7=1, and assume that A is small enough for the limits of integration to be extended to infinity. For 
the uncertainty A£, defined by (A£) 2 = f ((£\£ 2 \0) f ~ ( / ((CICIC» / ) 2 , we find 

A£ = \= . (67) 

We see that, for a fixed window width A, it is possible to make £ arbitrarily sharp when v is sufficiently large. On 
the other hand, when A is made small, in effect imposing the constraint of the clockless model, it becomes impossible 
to construct sharp eigenstates of £, reflecting the fact that, without the clock, there is no gauge-invariant operator 
corresponding to X or Given the coherent state \z))f , we can estimate the probability amplitude for finding the 
oscillator at the position £ at time r. Assuming that the limits of integration in ( |66| ) can be extended to infinity, 
which will be valid if v is large enough and if |£| < 1, where £ = JJz, we obtain 



'«£|*(r)))'~ exp 



(68) 



where ~ indicates the omission of a normalization factor and a phase. For a fixed window width A, this amplitude 
becomes sharply peaked at the classical trajectory in the semiclassical limit v — > 00. However, when A becomes 
very small, then (i) there is no longer a sharp peak at the classical trajectory and (ii) the probability becomes 
time-independent. Qualitatively, at least, we recover the situation in the absence of the clock, where there is no 
gauge-invariant operator corresponding to £ and all gauge-invariant amplitudes are time-independent (except possibly 
for trivial phase factors). 

It is apparent from these results that, as expected, we cannot fix the energy of the clock without losing both the 
observable £ and the time dependence which the clock was intended to reveal. However, it seems that one can find 
situations (for example, by taking v large, with IS.v fixed) in which both the energy of the clock and the eigenstates 
of £ are fairly sharply defined, and in which the notion of a test clock is therefore meaningful. 
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IV. DISCUSSION 



The absence of time-dependent observables in a time-reparametrization-invariant system can perhaps be understood 
by recognising that a complete description of a closed system is inevitably from the point of view of an observer external 
to the system. As viewed by such an observer, the system has, say, d + 1 dimensions, of which one corresponds to the 
"time" coordinate t (so the models considered in this paper have d = 0), and the "time" which passes inside the system 
is quite unrelated to any "time" experienced by this observer. For this observer, a state of the (d + l)-dimensional 
system corresponds to what an observer internal to the system might be supposed to regard as an entire history of 
a d-dimensional system. The observable quantities whose values distinguish one state of the system from another 
therefore characterise entire histories of the system, and cannot evolve with time in any straightforward sense. 

For a classical system which is not time-reparametrization invariant, this does not normally present a problem. 
Consider, for example, the model ([[]) with A set equal to 1. The general solutions to the equations of motion are of the 
form x(t) = X cos(wt) + (P/lu) sin(W) and q(t) = Q + lit. There are four independent observables, (X, P, Q, II), since 
a set of values for these quantities determines a history of the system. However, for a fixed value of t, the quantities 
x(t) and q(t), for example, are linear combinations of (X, P,Q,H) and are also observables. From such families of 
observables, parametrized by t, we can obviously construct functions such as x(t) which an internal observer might 
construe as representing the evolution with time of a single observable x belonging to his 0-dimensional "space" . This 
assumes that the internal observer, who does not appear explicitly in our model, does not significantly disturb the 
objects represented by x and q. If suitably equipped, he will be able to inspect two independent objects, an oscillator 
and a "free particle", associated respectively with the pairs of phase-space coordinates (X, P) and (Q,H). 

In the case of a time-reparametrizarion-invariant system, the families x(t) and q(t) have no gauge-invariant meaning, 
since they depend on t through the undetermined lapse function \{t). The gauge-invariant quantities which might 
qualify as observables are necessarily global quantities, in the sense that they cannot be associated with any particular 
value of t, and cannot evolve with time in any straightforward sense. For example, the gauge- invariant variables \ an d a 
in ( p^| ) and ( p0| ) specify the amplitude of the x oscillation and the phase of this oscillation corresponding to q = 0. One 
way out of this difficulty is to suppose that a specific lapse function is spontaneously selected, so that gauge invariance 
is broken, and gauge- variant quantities such as x(t) become physically meaningful. In general, this strategy seems 
unsound, since no physical principle gives rise to a preferred function X(t). In effect, the reparametrization invariance 
is removed artificially, and the constraint obtained by varying X(t) is then poorly justified. We note, however, that 
recent attempts to adapt the de Broglie-Bohm interpretation of quantum mechanics to quantum gravity ]T^ , p^ | do 
appear to produce a preferred foliation of spacetime and hence a natural definition of time evolution. 

It is tempting to account for time evolution in the following way. Consider the reparametrization-invariant model 
consisting just of a single oscillator (the model ([!]) with q omitted). The solution to its equation of motion is of the 
form 



x(t) = a cos (ip(t) +7) . (69) 

The amplitude a is constrained to have the value po, so there are no gauge-invariant variables, and the (0 + 1)- 
dimensional system has just one state available to it. Nevertheless, it seems intuitively clear that an observer inside 
the (0 + l)-dimensional universe would perceive an oscillation of the form X(r) = pqCOs(ujt), where r is the proper 
time along this observer's world line. Moreover, this assertion seems to be gauge- independent, since the lapse function 
serves only to relate r to an arbitrary coordinate t, which is irrelevant to the observer. Indeed, it is essentially this 
line of argument which enables an account to be given in elementary applications of general relativity of the time- 
dependent appearance of the universe to an observer. The discrepancy between the time evolution perceived by an 
internal observer and the single state apparent to an external observer arises from the fact that our model includes 
no description of the clock from which the internal observer gains his sense of time. 

To reconstruct what an internal observer's experience might be, it seems plausible that we must incorporate in 
our model at least a rudimentary description of this observer's clock. Classically, the variable q, whose behaviour is 
exhibited i n (120j ) mimics an instrument whose reading is linear in the geometrical proper time r. The quantity x(t) 
defined in (|2l[)defines a family of gauge-invariant observables which can be construed as giving the value of x which 
would be perceived by an internal observer when a proper time r has elapsed since his clock read zero. Moreover, 



the evolution of x(t) with r is governed by an equation of motion (25) of the standard Hamiltonian form. Quantum- 
mechanically, this equation of motion translates to a Heisenberg equation (|5^), again of the standard form. It thus 
seems consistent with the Copenhagen interpretation of quantum mechanics to suppose that an external observer 
might determine the state of the "universe" by measuring the value of x(ti), say with the result x\. One can then 
compute the probability | ( (a?2 3 T2|xi, t\))\ 2 that a measurement of xfe) will yield the value Xi- This probability may 
plausibly be interpreted as the probability that an internal observer, having determined the position of the oscillator 
as x\ will obtain the value xi from a measurement made after an interval T2 — T\ of his proper time. 
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To this extent, we recover precisely the usual formulation of time-dependent quantum mechanics. However, the 
status of the clock variable in this formulation requires further thought. In addition to x(t), we can construct from 
© the family of gauge-invariant variables q(r) — Qt, where Q — wpo sm X = ^VPo ~~ ~~ (-P/ 1 ^) 2 - Classically, 
the values of x(t) and q(r) can be determined simultaneously, and it might appear that these values represent the 
results of inspecting both the oscillator and the clock at time r. However, in contrast to the non-reparametrization- 
invariant system, the physical reduced phase space is now 2-dimensional. The two coordinates X and P are just 
sufficient to represent the position and momentum of the oscillator, and there are no further coordinates available to 
represent the clock as an independent dynamical object. As emphasised, for example, by Unruh Ej ], who has also 
considered the two-oscillator model, the variables a\ and 61, in ( p8|) and j2S|), which one would ordinarily want to 
treat as Schrodinger picture operators representing the oscillator and clock positions, do not correspond to operators 
on the physical Hilbert space. Indeed, the variable x{t) in (|2l])can be defined only when we solve the constraint, 
eliminating q as an independent dynamical variable. Correspondingly, the evolution of q(r) cannot be obtained from 
a Hamiltonian equation of motion, since Q commutes with Hq (equation (p6|)). Quantum-mechanically, x(t) and q(r) 
do not commute, so their values cannot be determined simultaneously. Moreover, the Hilbert space of this model is 
finite-dimensional, so the readings of any object we wish to regard as a clock must be drawn from a finite, discrete 
set of eigenvalues. From these considerations, it is apparent that a statement such as x(t±) = x\ cannot be taken 
as implying that both the oscillator and the clock have been inspected, yielding the value x\ for the position of the 
oscillator and the value T\ for the time. 

Thus, r cannot be regarded as a phenomenological time deduced by an observer from his inspection of a physical 
clock which reads q. Rather, the definition of x(t) reflects a decision on the part of an external observer (or of a 
theoretician) to study time evolution from a particular point of view. One might, for example, imagine an automatic 
observing apparatus, which performs a sequence of tasks under the control of an internal clock, q. To study time 
evolution "from the point of view" of this apparatus, it is natural to solve the Hamiltonian constraint by eliminating 
the unobserved quantity q, in order to arrive at the time-dependent observable x(r). In this sense, r represents a 
"Heraclitian" time of the kind sought by Unruh and Wald j|] , which "sets the conditions" for a measurement to be 
made. Like the parameter t in a non-reparametrization-invariant system, r itself is not a measurable quantity. 

Furthermore, it appears quite unnecessary for the eliminated variable to be linear, or even monotonic, in r. The 
physical Hilbert space associated with the model ( p7| ) and the algebra of the gauge-invariant operators X and P 
obtained by eliminating q is independent of the parameter N which determines the frequency of the q oscillator. The 
role of the clock q in passing from either ( |l9| ) and ( p0| ) or ( |32| ) and ( |33| ) to the gauge-invariant observable ( pl| ) is that 
it distinguishes a history in which, for example, x = pocosxeosa when q = from one in which some other value 
of q corresponds to this value of x. In essence, the definition of a gauge-invariant observable x(t) becomes possible 
when an origin r = can be specified in a coordinate-independent manner, and the role of q is to provide this origin. 

The view of time evolution described here is similar, though not identical, to that proposed by Rovelli In ||, 

Rovelli considers essentially the model (|27]) with N = 1. He defines a gauge-invariant classical observable given, in 
our notation, by 



x(s) = po cos x cos 



a + cos 1 



Po sm x 



(70) 



which is the value of x when q = s. Like our r, the clock time s can assume a continuous range of real values. (Note 
also that, like our x(r), this x(s) is gauge- variant at x = 0, since it depends on a, though in fact only x(s = 0) 
is well-defined.) The qualitative discussion given in || suggests that Rovelli wishes to regard x(s) as representing 
the evolution of x relative to the time recorded by a physical clock. For the reasons given above, however, it is not 
possible to regard s and x(s) as the two results of inspecting a clock and an oscillator simultaneously. 

The operator corresponding to ( |70| ) is approximately self-adjoint, and obeys an approximate Heisenberg equation 
of motion, only when restricted to a region of the physical Hilbert space which Rovelli calls the "Schrodinger regime" , 
corresponding roughly to states in which the wavefunction tp(x,q) is sharply peaked around the classical trajectory. 
For this reason, Rovelli contends that the notion of time does not exist at a fundamental level, but emerges only in 
a suitable semiclassical limit. By contrast, our operator x{t) is exactly self-adjoint, and obeys the exact Heisenberg 
equation (|53|), with respect to the variable r which classically corresponds to the geometrical proper time, and this 
provides a counterexample to Rovelli's contention. 

We would like, of course, to speculate that the view of time evolution proposed here can be extended to canonically 
quantized general relativity. In such an extension, it seems likely that the role of r would be played by the proper 
time along the trajectory of an observer for whose observations we wish to account, while x(r) would correspond 
to local observables defined on this trajectory. To define, say, time-dependent observables throughout a space-like 
hypersurface, one would presumably need to introduce a space-filling family of observers (or, at least, their clocks). 
These might correspond to a reference fluid of the kind described by Brown and Kuchaf |15| , though the relationship 
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of our interpretation to that proposed by these authors is not entirely clear to us. 
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Figure caption 

FIG. 1 The squared magnitude of the wavefunction (^) (circles) for a coherent state of the constrained oscillator, 
which is defined only at the discrete values Xj, compared with the corresponding Gaussian wavepacket (solid curve) 
for the standard unconstrained oscillator. The normalization of the Gaussian packet has been adjusted so that the 
peaks of both wavefunctions have the same height. Time evolution is depicted over a half cycle of oscillation, with 
amplitude slightly smaller than the largest eigenvalue Xj. In this case, v = 10. 
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